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ABSTRACT 

The emission spectra of TeV blazars extend up to tens of TeV and the emis- 



sion mechanism of the TeV 7-rays is explained by synchrotron self-Compton 
scattering in leptonic models. In these models the time variabilities of X-rays 



and TeV 7-rays are correlated. However, recent observations of 1ES 1959+650 
and Mrk 421 have found the "orphan" TeV 7-ray flares, i.e., TeV 7-ray flares 
without simultaneous X-ray flares. In this paper we propose a model for the "or- 
phan" TeV 7-ray flares, employing an inhomogeneous leptonic jet model. After 
a primary flare that accompanies flare-up both in X-rays and TeV 7-rays, radi- 
ation propagates in various directions in the comoving frame of the jet. When a 
dense region in the jet receives the radiation, X-rays are scattered by relativistic 
electrons/positrons to become TeV 7-rays. These 7-ray photons are observed as 
an "orphan" TeV 7-ray flare. The observed delay time between the primary and 
"orphan" flares is about two weeks and this is accounted for in our model for 
parameters such as F = 20, d = 4 x 10 17 cm, a = 3, and 77 = 1, where T is the 
bulk Lorentz factor of the jet, d is the distance between the central black hole 
and the primary flare site, a/T is the angle between the jet axis and the direction 
of the motion of the dense region that scatters incoming X-rays produced by the 
primary flare, and rj/T is the angle between the jet axis and the line of sight. 
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Subject headings: BL Lacertae objects: individual (1ES 1959+650, Mrk 421) - 
galaxies: active - radiation mechanisms: nonthermal 

1. Introduction 

Blazars are a subclass of active galactic nuclei and their high energy emission is generated 
in regions near the central super massive black holes. Because of their intense and rapidly 
variable radiation, the emission regions are thought to be in relativistic jets closely aligned 
to the line of sight (see Aharonian 2004; Krawczynski 2004, for a review). The spectral 
energy distributions (SEDs) of blazars are characterized by double broad peaks in the v-vF v 
representation. One peak is in the optical to X-ray regions and the other is in the 7-ray 
region. Very high energy (100 MeV - 1 GeV) 7-rays from blazars were discovered by EGRET 
(Energetic Gamma- Ray Experiment Telescope) on board Compton Gamma-Ray Observatory 
(Hartman et al. 1999). Recent observations by Cerenkov telescopes have revealed that many 
blazars emit 7-rays up to tens of TeV (Aharonian et al. 2005, for a review). Flare activities 
in X-rays and 7-rays from blazars are also known. Although the emission mechanisms of 
blazars are still under study, the SEDs of blazars are well explained by leptonic or hadronic 
models. In the framework of the leptonic models, the emission in the optical to X-ray region 
is attributed to synchrotron radiation by nonthermal electrons/positrons in the jet and the 
TeV 7-rays are accounted for by the synchrotron self-Compton (SSC) model (Maraschi et al. 
1992). The SSC model assumes that synchrotron photons are inverse- Compton scattered 
by the same nonthermal electrons/positrons that emit synchrotron photons. If the SSC 
model is applied, the time variabilities of the X-rays and 7-rays should be correlated (e.g., 
Mastichiadis & Kirk 1997; Li & Kusunose 2000; Kusunose et al. 2000), and indeed most 
flares occur almost simultaneously in X-rays and 7-rays. However, the recent observations 
of 1ES 1959+650 (Krawczynski et al. 2004; Daniel et al. 2005) found that there is a TeV 
7-ray flare that is not accompanied by a X-ray flare, which was observed on June 4, 2002 
(MJD 52,429) by the Whipple telescope. This is called the "orphan" TeV 7-ray (OTG) 
flare. The OTG flare occurred about 15 days after a usual flare in which X-ray and 7-ray 
flares were observed contemporaneously. More recently another OTG flare might have been 
observed from Mrk 421 (Blazejowski et al. 2005), although the X-ray flux seems to have 
peaked about 1.5 days before this TeV flare and this may not be a true OTG flare. 

The existence of the OTG flares is challenging for the leptonic models. Bottcher (2005) 
has recently proposed a model for the OTG flare of ES 1959+650, i.e., the hadronic syn- 
chrotron mirror model. He assumed that a fraction of X-rays produced in the primary flare by 
synchrotron radiation of leptons is reflected by a plasma cloud with scattering depth ~ 0.1 lo- 
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cated in the direction of the jet propagation. The reflected X-rays collide with protons in the 
jet and pions are produced. Subsequently a large number of neutral pions decay into 7-rays, 
which are expected to be observed as TeV 7-rays. If this model is applied to 1ES 1959+650, 
it is found that the reflecting plasma is located at a distance ~ 3(r/10) 2 At 2 o P c from the 
central black hole, where T is the bulk Lorentz factor of the jet and At2o — At/(20days) 
is the normalized time delay between the primary and OTG flares. This model requires 
unreasonably large values of proton density in the jet and the hadronic jet power (Bottcher 
2006: erratum). Thus the hadronic synchrotron mirror model may not be applicable to 
OTG flares, although it may be still viable if the effects that reflected synchrotron photons 
increase as the blob approaches to the mirror are taken into account. (More recently Reimer 
et al. (2005) calculated the predicted neutrino spectrum due to the decay of charged pions 
in the framework of the hadronic synchrotron mirror model.) 

In this paper, we propose another model of the OTG flares in the framework of the SSC 
model, assuming that the emission region is not homogeneous, which is different from the 
conventional SSC model. We assume that the injection of nonthermal electrons/positrons 
triggers the primary flare where both X-rays and TeV 7-rays are emitted. If the jet is not 
uniform but consists of a few patchy regions, X-rays that were produced in the primary flare 
impinge on another dense region of the jet. This sudden increase of X-ray photons results 
in strong TeV 7-ray flux by inverse Compton scattering, which is expected to be observed 
as an OTG flare. Because there is a time delay between the primary flare in one region and 
the scattering in another region, TeV 7-rays are observed as an OTG flare. Note that TeV 
7-rays from the primary flare are not scattered in the OTG flare site because of Klein-Nishina 
effects. The observed delay time between the OTG flare and the previous flare is about two 
weeks in 1ES 1959+650 (Krawczynski et al. 2004). The delay time between the primary and 
OTG flares in the comoving frame of the jet is much shorter than the synchrotron cooling 
time of nonthermal electrons/positrons (see §3 for detail), if the magnetic field in the jet is 
about 0.1 G. Although high energy electrons/positrons injected into the primary flare are 
rapidly cooled, nonthermal electrons/positrons that emit TeV 7-rays in a quiescent state may 
be still continuously injected into the jet and they may contribute to the secondary flare. 
In this paper we do not specify the acceleration mechanisms of those electrons/positrons or 
calculate the emission spectrum of the flare but focus our work on the kinematics of the jet. 

The kinematics of the jet is described in §2 and the jet properties are presented in §3. 
A summary of our results and discussion are given in §4. 



-4- 



2. Kinematics of Jets 

We first study the kinematics of a jet shown in Figure 1 (Case 1). We assume that a 
jet (a patchy plasma shell) with bulk Lorentz factor T is ejected radially from point O that 
is close to the central black hole. The observer is at angle # b s = t}/Y from the jet axis, 
where 77 ~ 0(1). At P the primary flare occurs, which is triggered by, e.g., the internal 
shock. Here we assume that the high energy particles are confined in a patchy shell and 
high density regions are located at P and Sp in the shell. The emission mechanisms in a 
flare are synchrotron emission and SSC for X-rays and TeV 7-rays, respectively. A fraction 
of the emitted photons propagate along directions that are slightly off from the jet axis and 
they are injected into the high density region around S, where S was located at Sp when 
the primary flare occurred. The injected TeV 7-rays are in Klein-Nishina regime and are 
not scattered but X-rays are Compton-scattered to become TeV 7-rays. This triggers the 
secondary (OTG) flare in the TeV region. 

We assume that the photons emitted at P reach S after time interval 5t in the rest frame 
of the black hole and that the propagation direction of those photons is denoted by 8f. For 
the geometry shown in Figure 1, equations to be solved for St and Of for given d, £, and T 
are as follows: 

(d + (3c5t) sin£ = c5t sm6 f , (1) 

and 

(d + (3c5t) 2 = d 2 + (c5t) 2 - 2 dc5t cos(n - 9 f ), (2) 

where c is the light speed and (3 = a/1 — 1/T 2 . Here we assumed that the distances between 
O and P, OP, and O and Sp, OSp, are the same, i.e., OP = OSp = d. Next we assume that 
£ = a/T, where a ~ Oil) is a constant. Using equations (1) and (2), we obtain the equation 
for x = cSt/d: 

x 2 + 4TVf^T f—-x)+ 1^5! - Aa 2 + 4(a 2 - l)r 2 = 0. (3) 
\ x J x 2 

There are four solutions for x, but the physically allowed solution is unique and given by 
x = rVf 2 ^! - Vr 4 - (1 + a 2 )T 2 + a 2 

4 2r - a(r - 1)+ r-(UW' +a T/- (a - r) } ■ (4) 

When r ^> 1, c5t/d is approximated as 

~d~^Y + V a +T + G( } - (5) 
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From equation (1) we obtain 

/ r 



a + 



+ V4 + a 2 



CK 



(6) 



Because the angle between the jet axis and the line of sight is # b s , the observed delay time 
of the secondary flare from the primary flare is given by 

5t ohs = [1 - 008(0/ - 9 ohs )] 5t^^ f(a) d - , (7) 

where 



1 / 2 / „ . a 4 



2 + a 2 + V4a 2 + a 4 - an 4 + + 2n 2 

a 2 + ^Aa 2 + a±; 



(8) 

The behavior of f(a) is shown in Figure 2. For example, / ps 38.2 for a = 3 and rj — 1. 
Note the strong dependence of / on a. The numerical value of <5t b s is estimated as 



20\V d 



fobs ~ 15 days is obtained for V = 20, d — 4 x 10 17 cm, 77 = 1 , and a = 3. This value 
is close to the observed delay time of the OTG flare in 1ES 1959+650 (Krawczynski et al. 
2004). These parameters give the distance of S from the central black hole d + f3c8t ps 1.4 
pc. 

There is another possibility of the jet geometry (Case 2) as shown in Figure 3. Here 
we assume that the X-ray source is located at P that is away from the jet axis and that the 
secondary flare occurs at Sp on the jet axis. The solution for c5t/d is the same as equation 
(5) and 9f is given by equation (6). If a flare at Sp occurs simultaneously with a flare at P 
and the flare at Sp is observed, the observed time delay between the flares at Sp and S is 
given by 

5t bs = (1 - /3 cos 9 ohs )5t ~ =2/1(0;) -, (10) 

1 z c 

where 



The values of f\(a) for rj — 0, 1, and 2 are shown in Figure 4. If T = 20, a = 3, rj = 1, 
and d = 4 x 10 17 cm are assumed, the delay time is about 4 days and too short to explain 
the observed delay time. If 77 is larger, the delay time becomes larger, but the observed 
luminosity decreases significantly because of the beaming effects. 
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If the primary flare at P is observed and a flare at Sp is not observed, the delay time 
between the flares at P and S is given by 



[1 - cos(£ obs - f + e f )]8t » X f 2 (a) | if f < 6 obs 
8t obs = { 1 ^ 1 (12) 

' ;i- cos e f )8t « ^ fc(a) § ife>^obs, 



where 

/ 2 (a) = 4m? + a 2 (r? 2 - 1) + oV4 + a 2 (l + ?? 2 ), (13) 

and 



a 4 + 3a 2 , 1 + a 2 / 2 , a 4 
/s(a) = ^ + _ 2~ V T" ^ 

/ 2 (a) for 77 = 1.5 and 2 and /3(a) are shown in Figure 5. When £ < 6> obs , ?7 > 1.5 and a > 1.5 
are needed to obtain 5t obs ~ 2 weeks, if T = 20 and d — 4 x 10 17 cm are assumed. When 
£ > #obs, « > 2.7 is required to explain 5t ohs ~ 2 weeks for the same parameter values as the 
above. 



3. Properties of the Secondary Flare 
3.1. Timescales 

The solution obtained in §2 gives the distance, SpS, that the emission region travels 
during St as 

'a 2 H a* 



f3 cS t~ [— + ^ + -) d . (15) 



In the comoving frame of the jet, the travel time is given by 



If a = 3, we obtain 8t' ~ 574 days. Here primed quantities are evaluated in the comoving 
frame of the jet. 

Synchrotron cooling rate per electron in the comoving frame of the jet is given by 
P s ' yn = (4/3)o"tc/5 /2 7 /2 -u^, where u' B is the magnetic energy density. Then we obtain the 
cooling time as t' syn = 7' 'm e c 2 / 'P s ' yn . The characteristic energy of synchrotron radiation is 
given by e' syQ = SheB'^' 2 / (2m e c) and the observed synchrotron photon energy is Doppler 
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shifted, i.e., e S yn ~ r6g yn . Because synchrotron photons with energy c' S y n are emitted by 
electrons with a Lorentz factor 

we obtain 

^~ L7 (o^)" 3/2 S) 1/2 (T5Sv)" 1/2 d ^ (18) 

Then we find that t' syQ < St'. Because the cooling by inverse Compton scattering is not 
negligible, the cooling time is much shorter than t' syn . Thus relativistic electrons/positrons 
injected into the primary flare are already cooled when the secondary flare occurs. For the 
secondary flare to be caused, the injection of relativistic electrons/positrons is necessary. 
The injection rate, however, does not need to be as high as in the primary flare. When 
synchrotron photons from the primary flare impinge on the dense region, which is located at 
the right place as described in §2, the secondary flare occurs by inverse Compton scattering 
of those photons. 



3.2. Beaming effects 

There is a slight velocity difference between Sp and P. The relative velocity between 
points Sp and P normalized by the light speed is given by f3 Te \ = 2(3 sin(£/2) ~ a[3/T. Thus 
the relative velocity is negligible for T ^> 1 and points Sp and P are safely assumed to be 
in the same inertial frame. The observed luminosity of TeV 7-rays from S is enhanced by a 
factor D| because of the relativistic effects in Case 1, where 

Di = r[l-/3cos(£-0 o bs)] (19) 
2 r , ( a -r,)* + 6(a-r,) 2 -3 1 , ^_ 3 , 



~ i + (a-vr T+ 6[i + (Q -,)T r + 0(r (20) 

In Case 1, the ratio of the luminosity enhancements by the beaming effects at S and P is 
r = (Zy/W 4 = {(1 + V 2 )/[l + (« - where 

n 1 2 ^ + 6^-3 1 3 

r(l - f3cos6 ohs ) 1+r] 2 6(1 + r] 2 ) 2 T 
If a = 3 and 77 = 1, r w 1/40. In Case 2, r = (D ohs /D ( ) 4 w 40. 
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4. Summary and Discussion 

In this paper we have proposed a structured jet model for the OTG flares. The high 
energy emission spectra from TeV blazars have been well explained by the leptonic jet models. 
Because the leptonic jet models assume synchrotron self-Compton scattering for the emission 
mechanism of very high energy 7-rays, it is expected that flares occur in both X-ray and 
7-ray energy bands almost simultaneously. However, the OTG flares were not accompanied 
by the simultaneous increase in the X-ray flux (Krawczynski et al. 2004; Blazejowski et 
al. 2005), and this is thought to be a challenging issue for the leptonic jet models. Our 
model of a structured jet assumes that 7-rays are emitted in different components in the 
jet. X-rays produced by a flare in a region, where flares in X-rays and TeV 7-rays occur 
simultaneously, are scattered in another region of the jet and the scattered photons are 
observed as a secondary (OTG) flare. Note that TeV 7-rays of the primary flare are not 
scattered by electrons/positrons because of Klein-Nishina effects. The emitted 7-ray energy 
is estimated as e' ssc = 7'm e c 2 , because the scattering occurs mainly in Klein-Nishina regime. 
Then the observed 7-ray energy is roughly given by 

fssc ~ r 4sc ~ 5.5 x 10" g) - 1 ' 2 (I) V2 ( Jg^) " 2 eV, (22) 
using equation (17). 

Because the X-rays emitted in the primary flare need a propagation time before they are 
scattered in different parts of the jet, this results in the time delay between the primary and 
secondary flares. In our model, the delay time between the primary flare and the secondary 
TeV 7-ray flare is dependent on T, a, 77, and d. The observed delay time is about 15 days 
for 1ES 1959+650 and this value is obtained for V = 20, a — 3, 77 = 1, and d = 4 x 10 17 
cm in Case 1. As the value of a decreases the value of d increases for a given value of <5t bs- 
If the value of T is known from the observations of the primary flare, the values of a and d 
might be reasonably guessed from the delay time. It is not clear from the numerical values 
estimated in this paper whether Case 1 is favored to Case 2 or vice versa. Rare observations 
of OTG flares might be due to that this phenomenon needs a specific structure of the jets 
as shown in Figures 1 or 3. 

The duration of the secondary flare depends on the time interval of the primary flare 
and the angular size of the secondary flare site, A9. For example, in Case 1, if A6 is given by 
Aa/T and the primary flare is impulsive, the observed duration of the second flare is given 
by the difference between the values of <5t b s for a and a + Aa: 

/20\V d \df(a) A , , , 
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where Aa <C 1 is assumed. The values of g(a) = df /da are shown in Figure 6 for Case 1. 
For example, g m 54.4 for r] — 1 and a — 3. If the observed duration of the secondary flare 
is about 5 hours, Aa ~ 0.01 for a = 3, T = 20, and d — 4 x 10 17 cm. The observed duration 
of the secondary flare may depend on not only the value of Aa but also the duration of the 
primary flare, and the effect of the latter elongates the duration of the secondary flare. 

We have shown that the time lag between the primary and secondary flares are well 
accounted for by our model. Next we argue that the density of soft photons (X-rays) from 
the primary flare is high enough to produce the OTG flare if the particle injection time in the 
primary flare is sufficiently long, and that the particle injection rate needed in the secondary 
flare is not unreasonably high. 

First, it should be noted that photons from the primary flare blob (PB) are not injected 
from the rear side of the secondary flare blob (SB). In the comoving frame of the blobs, the 
angle between the line of sight and the incident direction of photons from PB is about 90 
degrees. Then the energy boost by inverse Compton scattering in SB should be effective, 
although the detailed calculations of radiative transfer are necessary to obtain the emission 
spectrum of TeV 7-rays from SB. 

Next, we discuss the amount of radiation received by SB and the cooling time of leptons. 
In the following, primed symbols are quantities during the primary flare in the comoving 
frame and doubly primed symbols are quantities during the secondary flare in the comoving 
frame. We assume that PB and SB are spherical plasma clouds for simplicity. The radii 
of PB and SB in the comoving frame at the occurrence of the primary flare are denoted 
by -Rp rim and R' sec , respectively. The radius of SB at the time when the secondary flare 
occurs is denoted by R'^ c - The solid angle extended by SB for a point in PB is given by 
^scc ~ ^-^sec/^sec where d' scc is the distance between PB and SB: d' scc ~ £d = (a/T)d. 
If SB expands proportionally with the distance from the central black hole, Q' sec = ^ s ' ec . 
If radiation from PB is emitted isotropically in the comoving frame, the fraction of the 
radiation received by SB is given by 

__n^ = lT^RZ 1 /r\ 2 /3\ 2 / i4 c \ 2 / 4xl0 17 cm \ 2 
rrad ~ 4tt 4 a 2 d 2 ~36V20; \a ) V2xl0 16 cmJ V d J ' 1 ' 

We assumed R' scc ~ 2 x 10 16 cm; this value is close to the values assumed for the emission 
region of the flare models of 1ES 1959+650 in Krawczynski et al. (2004). Thus the soft 
photon (X-ray) energy density in SB is roughly about 1/30 of that in PB and the Comp- 
ton cooling time of nonthermal particles in SB is longer by 30 times than in PB. On the 
other hand, the flare duration in PB is determined by the injection duration of nonthermal 
particles, because the cooling time should be short owing to the large energy densities of 
soft photons and magnetic fields. Thus the observed OTG flare is explained, if the particle 
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injection duration in PB is about 30 times longer than the cooling time in SB. Furthermore, 
the energy densities of magnetic fields and internal synchrotron photons in SB are smaller 
than in PB, so that emission other than TeV 7-rays is too weak to be observed. 

In the following we estimate the required injection rate of nonthermal leptons in the 
secondary flare. The observed energy flux of the secondary flare at 600GeV is vF v ~ 3 x 10~ 10 
ergs s^ 1 cm" 2 . The luminosity distance to 1ES 1959+650 (z = 0.047) is d L = 210 Mpc, 
assuming that H = 70 km s _1 Mpc -1 , — 0.7, and fl m = 0.3, where H , and Q m are 
the Hubble constant and the density parameters of dark energy and matter, respectively. 
Then the luminosity of the very high energy gamma-rays is Lvhe.scc ~ 1.6 x 10 45 ergs s _1 . 
The luminosity in the comoving frame is given by 

^VHE,sec „ „ ln 39 / ^scc \ _i 



L„ v n£j,scc „ „ in39 / ^ sec \ — i / nr \ 

VHE,sec = —jyT 9 ' 9 X 10 J er S SS > ( 25 ) 

where X> sec is the beaming factor of SB and V scc = 8 and 20 for Cases 1 and 2, respectively, 
if T = 20, a — 3, and 77 = 1 . The distance between the central black hole and SB is 
d sec = d + (3c5t ~ lid, where a = 3 is assumed. Assume that SB expands linearly with the 
distance from the central black hole. Then the radius of SB is given by 

R'^ = %£ R' scc ~ Hi4 c = 2.2 x 10 17 ( - ) cm. (26) 

d V 2x10 cm / 

We obtain the energy injection rate per unit volume during the secondary flare: 

C ^ W ^ 22 x l °- 13 (2^) " S (If) " -" 3 < 27 > 

The injection rate of nonthermal leptons in SB is also given by 

N „ _ L VHE, S ec _ 9 4 1Q 40 ( ^scc\ ( (t) \ -1 ^gx 

where (7) is the value of the characteristic Lorentz factor of injected leptons in SB. When 
the energy injection rate during the primary flare is denoted by q' prim , we obtain 

<?s'ec -^VHE.scc/ (if^sec) [^^]vHE,sec ( ^prim \ ^ ( -^prirn \ 3 / on N 



9prim -^VHE.prim/ (lf-^>rim) VHE,prim \ Ace / \ R'L 

where [^i^]vHE,prim and [zaF„] vHE,sec are the observed quantities of vF v in the primary and 
secondary flares, respectively. If i?g CC ~ -R pr i m and R" ec ~ H-Rgec are assumed, q" ec /q' w - nn ~ 
Ci [zaF^] V HE,sec/[^^]vHE,prim, where C\ ~ 3.0 x 10~ 2 and 1.9 x 10~ 5 for Cases 1 and 2, respec- 
tively. From the observed data (Krawczynski et al. 2004), we may set 
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[^]vHE,sec/[^]vHE, P rim ~ 1, and g" ec ~ cigp rim is obtained. Thus the required value of q" ec 
is not unreasonably large. This particle energy is transfered to X-rays by inverse Compton 
scattering to emit TeV 7-rays in the secondary flare. 

Our model assumes a structured jet and the occurrence of OTG flares depends on 
the geometry of dense regions in the jet. Alternatively the dense regions might be widely 
distributed, e.g., the layer-spine model by Ghisellini et al. (2005), and the value of a depends 
on where the region of particle acceleration at pc-scale exists. 

Finally we comment on the observed large scale structure of jets. The line of sight is 
inside the opening angle of jets in our model, and this may result in a core-halo structure 
of the emission region. However, since the emission regions of high energy radiation are 
located very close to the central black hole, the large scale structure of jets is affected by 
the environment of the jets. That is, the interaction between jets and gases surrounding 
the central region may cause the bending of the jets. It is also possible that the opening 
angle may change with the distance from the central region. These may account for the jet 
structure observed by VLBI. 

We thank the anonymous referee for the valuable comments. This work has been par- 
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Ministry of Education, Culture, Sports, Science and Technology of Japan. 
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Fig. 1. — Case 1. (a) The jet is ejected from O and consists of patchy regions, and they 
move radially with Lorentz factor T. The primary flare occurs at P. During time interval 5t, 
the dense region around Sp moves to S and a fraction of X-rays emitted at P travels toward 
S. Then the secondary ("orphan") flare is triggered at S. The observer is located at angle 
#obs = 77/r off the jet axis, and the points S P and S are on the line at angle £ = a/V off the 
jet axis. The distances OP and OS P are d, S P S = (3c5t, and PS = cSt. (b) Emission regions 
in the comoving frame of the jet. X-rays emitted in the primary flare are injected into S and 
they are scattered by nonthermal particles in S. 
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Fig. 2. — f(a) for r) = (solid), 0.5 (dashed), and 1 (dot-dashed). 
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Fig. 3. — Case 2. Same as Fig. 1 except that X-rays come from the off-axis region around 
P. 
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Fig. 4. — fi(a) for r\ = (solid), 1 (dashed), and 2 (dot-dashed). 
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Fig. 5. — /2(a) for r] = 1.5 (dot-dashed) and 2 (dot-dot-dashed) and /3(a) (solid). /2(a) 
and /3(a) are used to calculate 5t ^ s in Eq. (12) for a < rj and a > r], respectively. 
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Fig. 6. — g(a) for r\ = (solid), 0.5 (dashed), and 1 (dot-dashed). 



